Quantum Melting of Spin Ice: Emergent Cooperative Quadrupole and Chirality 
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A quantum melting of the spin ice is proposed for pyrochlore-lattice magnets PT2TM2O7 (TM =Ir, 
Zr, and Sn). The quantum superexchange Hamiltonian having a nontrivial magnetic anisotropy 
is derived in the basis of atomic non-Kramers magnetic doublets. The ground states exhibit a 
cooperative ferroquadrupole and pseudospin chirality, forming a magnetic analog of smectic liquid 
crystals. Our theory accounts for dynamic spin-ice behaviors experimentally observed in Pr2TM207. 

PACS numbers: Valid PACS appear here 



It has been a great challenge to realize unconventional 
spin-liquid states in three-dimensional magnets. It is 
achieved by preventing a dipole long-range order (LRO) 
of magnetic moments, which requires appreciable quan- 
tum spin fluctuations and geometrical frustration of mag- 
netic interaction The importance of the geometri- 
cal frustration is manifest as in pyrochlore systems [IrlH . 
In particular, in the dipolar spin ice ii^T^Oy (i?=Dy 
or Ho) [fHH , the rare-earth magnetic moment located at 
each vertex of tetrahedrons points either inwards ( "in" ) 
to or outwards ("out") from the center (Fig. Q] (a)). 
The nearest-neighbor ferromagnetic coupling mainly due 
to the magnetic dipolar interaction favors macroscopi- 
cally degenerate "2-in, 2-out" configurations without any 
LRO, forming a magnetic analog of the water ice [8|. 
Then, the classical spins are quenched into one of the de- 
generate ground states Usually, a thermal heating is 
required for melting the quenched spin ice. Here, we pur- 
sue an alternative possibility that quantum fluctuations 
melt the spin ice: the quantum entanglement among de- 
generate states lift the macroscopic degeneracy, suppress 
the spin-ice freezing, and lead to a distinct ground state. 

A realistic approach to the quantum melting of the spin 
ice is to choose a rare-earth ion with fewer / electrons 
and a smaller magnetic moment, e.g., Pr 3+ . In rare-eath 
ions with fewer f electrons, the 4/ wavefunction is spa- 
tially extended [lOj and can then be largely overlapped 
with the O 2p orbitals at the 01 site (Fig. Q] (a)) in the 
pyrochlore lattice. Besides, for Pr 3+ ions, the magnetic 
dipolar interaction, which is proportional to the square of 
the moment size, is reduced by an order of magnitude to 
0.1 K between the nearest-neighbor sites, in comparison 
to that for Dy 3+ ions. Then, the superexchange interac- 
tion due to virtual f-p electron transfers, which provides 
a source of the quantum nature, is expected to play cru- 
cial roles in P^TA^Oy (TM: a transition metal). 

Recent experiments on P^S^Oy [HI], P^Z^Oy [13], 
and Pr2lr207 [l3( have shown that the Pr 3+ ion provides 
the (111) Ising moment described by a non-Kramers mag- 
netic doublet. As_in_the spin ice, any magnetic dipole 
LRO is absent 



FIG. 1: (Color) (a) Pr 3+ ions (red) form tetrahedrons (dashed 
lines) centered at O 2- ions (01) (blue), and are surrounded 
by O 2- ions (02) (blue) in the D^d symmetry as well as by 
TM ions (green). Each Pr magnetic moment (bold arrow) 
points to either of the two neighboring Ol sites, (b) The 
local coordinate frame (x r ,y r , z T ) from the top. Upward and 
downward triangles of the O 2- ions (02) are located above 
and below the hexagon of the TM ions, (c) The Pr pyrochlore 
lattice. The phase <p r , r i in Eq. ([2]) depends on the color of the 
bonds. The global coordinate frame (X, Y, Z) is also shown. 
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16]. Pr2lr207 shows a metamagnetic 
transition only when the magnetic field is applied in the 
[111] direction [13], indicating the ice-rule formation due 



to a ferromagnetic coupling J ~ 1.4 K [1J]. On the other 
hand, the Curie- Weiss temperature Tew is antiferromag- 
netic for the zirconate 12] and iridate [l3j . unlike the 
spin ice. The stannate shows a significant level of low- 
energy short-range spin dynamics [la ], which is absent in 
the classical spin ice. Furthermore, the iridate shows the 
Hall effect at zero magnetic field without magnetic dipole 
order [3] , suggesting an onset of a chiral spin liquid Q 
at a temperature ~ J due to quantum fluctuations. 

In this Letter, we derive the realistic effective model for 
Pr 4/ moments on the pyrochlore lattice. It contains ap- 
preciable quantum nature leading to a cooperative ferro- 
quadrupolar ground state, accompanied by crystal sym- 
metry lowering from cubic to tetragonal and a frustration 
in the chirality ordering. Our scenario explains unusual 
magnetic properties observed in Pr 2 TM 2 07 suggesting a 
dynamically fluctuating spin ice [l4l - [lq] . 

We start with f 2 configurations for Pr 3+ forming the 
tetrahedron centered at the 2 ~ ion (01) in Pr2TAi207 
(Fig. [1]) . The LS coupling gives the ground-state mani- 
fold 3 #4. Each Pr 3+ ion is placed in a crystalline electric 
field (CEF) which has the D 3 d symmetry about the (111) 
direction toward the 01 site. It is useful to define the lo- 
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FIG. 2: (Color online) (a) Local level scheme for / and p 
electrons, and the local quantization axes z r and z r i . (b) f-p 
transfer integrals, (c) Virtual hopping processes, n (n') and 
£ in the state f n p e f n represent the number of / electrons at 
the Pr site r (r r ) and that of p electrons at the 01 site. 



cal quantization 0jX1S Zj> 3jS this direction, as well as x and 
y axes as x r and y r depicted in Figs. Q] (a) and (b). The 
CEF favors J z = ±4 configurations for the total angu- 
lar momentum, which are linearly coupled to J z = ±1 
and =f2 because of the D$d CEF 18]. This leads to the 
atomic non-Kramers magnetic ground doublet 

\a z ) = a\J z = Aa z )+f3a z \J z = <J z )-j\J z = -2a z ), (1) 

with small real coefficients (3 and 7 as well as a = 
\Jl — (3 2 — 7 2 . The pseudospin a z = ± represents the 
direction of the Ising ("in" or "out") magnetic dipole 
moment, in contrast to the case of a nonmagnetic dou- 
blet labeled by the atomic quadrupole moment I19| in 
materials having other CEF symmetries, PrFe4Pi2 [20j, 
UPt 3 [21], and UPd 2 Al 3 For Pr 2 Ir 2 7 , the first ex- 
cited crystal-field level is a singlet located at 168 K and 
the second is a doublet at 648 K [l8[ . They are similarly 
large for Pr 2 Sn 2 07 [l5[. These energy scales are two or- 
ders of magnitude larger than our relevant energy scale 
J ~ 1.4 K. Hence we neglect these CEF excitations. 

Now we derive the effective Hamiltonian through the 
fourth-order strong-coupling perturbation theory. Vir- 
tual local f 1 and / 3 states have an energy gain of the 
Coulomb repulsion U and cost of 2U, respectively (Fig. [5] 
(a)), where the LS coupling has been ignored in com- 
parison with U . Creating a virtual p hole decreases the 
energy by the p electron level A compared with the f 1 
level. The f-p electron transfer is allowed only within the 
orbital l z — and ±1 manifolds, whose amplitudes are 
given by Slater-Koster parameters V p f a and V p f n [23| . 
respectively (Fig. [5] (b)). The local coordinate frames for 
the nearest-neighbor Pr sites are crucially different; for 
instance, z r ■ z r > = — 1/3. The perturbation expansion in 
Vpf a and VpfK is then carried out by taking into account 
the different local coordinate frames and the virtual pro- 
cesses (Fig. [5] (c)). The projection of this superexchange 
Hamiltonian onto the subspace of doublets (Eq. ([1])) leads 



to the pseudospin- 1/2 Hamiltonian; 

n.n 

H eS = J ^2 [<JrO- z , +25 (a+o-- +0--0-+) 

(r,r') 

+2q(e i ^'a+a+,+h.c.)], (2) 

wither* = (of±ier^)/2 and (ct^,ct^,ct^) = a r -(x r ,y r , z r ), 
where a r represents the Pauli matrix for the pseudospin 
at the site r. We have adopted the simplest gauge where 
the phase ip r y takes 0, 2ir/3, —2ir/3 depending on the 
color of the bond directions shown in Fig.Q](c), by rotat- 
ing the x and y axes by 7r/12 from those shown in Figs.[T] 
(a) and (b). This phase can not be fully gauged away, 
because of the noncollinearity of the (111) magnetic mo- 
ments and the three- fold rotational invariance of (r,a r ) 
about the [111] axes. Only a z contributes to the magnetic 
dipole moment J z , while a% ,v the atomic quadrupole mo- 
ment J Z J*' V , as can be shown by direct calculations. For 
a realistic case —0.37 < Vpf^/Vpfu < —0.02, the Ising 
coupling J between the nearest-neighbor pseudospins is 
found to be positive, i.e., antiferroic. This indicates the 
"ferromagnetic" coupling between the nearest-neighbor 
4/ magnetic moments because of the tilting of the local 
z axes, z r ■ z r i = —1/3. Then, it can provide a source of 
the ice-rule formation. The D 3 d CEF creates two addi- 
tional quantum-mechanical interactions; the pseudospin- 
exchange and pseudospin-nonconserving terms. Their 
coupling constants S and q are insensitive to U /V p f a and 
A/Vpf a but strongly depends on (3 and 7. Figure [3] (a) 
shows S and q as functions of /3 for the trigonal CEF 
with by keeping the ratio 7//? = 3. Henceforth, we 
adopt rough estimates U/V p f a = 5, A/V p f a = 4, and 
Vpfir/Vpfa- = —0.3 from first-principles calculations (pub- 
lished elsewhere), and (3 = 7.5% and 7 = 3/3 from the 
CEF analysis based on inelastic neutron-scattering ex- 
periments 18]. Then, we obtain S ~ 0.51 and q ~ 0.89, 
indicating the appreciable quantum nature. The two cou- 
plings play crucial roles in inducing a cooperative ferro- 
quadrupolar order instead of the classical spin ice [1] or 
the U(l) spin liquid [H]. 

A mean-field analysis @ on Eq. ([2]) yields two dis- 
tinct states, (i) Magnetic dipolar states characterized by 
a nonzero (o~ z ) have the lowest energy — 23 per tetrahe- 
dron at the wavevector q — ^-(hhl) with a being the side 
length of the unit cube (Fig. Q] (c)). (ii) A quadrupolar 
state with a nonzero (o~r' v ) has the energy — 2(5 + 2q) J at 
q = for 5, q > 0. Thus, for S + 2q > 1 as in our case, the 
atomic quadrupole moment o~% ,v should form the LRO 
without any dipole LRO. However, we will show below 
that the ground state is further replaced with a cooper- 
ative ferroquadrupolar state because of the quantum in- 
terplay between atomic dipoles a z and quadrupoles a^. ,v . 

First let us solve Eq. <j2j) on a single tetrahedron. The 
similar analysis on a distinct model for Tb 2 Ti 2 C>7 [25| 
has been employed to discuss the RVB-singlet quan- 
tum spin ice [26j . With increasing [3 and thus 7 from 
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FIG. 3: (Color online) (a) Coupling constants 8 and q versus 
/3(= 7/3). The arrow points to the experimentally estimated 
value of /3 [181 ]. (b) Outward normal vectors (green arrows) 
of the surfaces of the tetrahedron, used to define the chirality 
kt- (c) Solid angle subtended by four pseudospins a ri . (d) 
Distribution of the tetrahedral magnetic moment Mt in the 
cooperative ferroquadrupolar state ({Qt Z ) > 0). The arrows 
represent the lattice deformation linearly coupled to Q zz ■ 



0, three classical levels corresponding to the "2-in, 2- 
out", "3-in, l-out"/"l-in, 3-out" (AE = 2 J), and "4- 
in" / "4-out" (AE = 8 J) configurations are split to three 
doublets, three triplets, and one singlet. In our case, 
the ground-state manifold has the E g symmetry with 
the double degeneracy x = ± as described as = 
71 Er=±( el¥x l^0 +e- l ^ x \rY) + \tZ)) + c 4 |4 X ) with 
real coefficients ci and oj. Here, the orthonormal state 
| + 4)/| — 4) represents the "4-in" /"4-out" configura- 
tion, while | ± X), | ± Y), and | ± Z) denote the 
"2-in, 2-out" having the net magnetic dipole moment 
Mt = M J2r a r^r, pointing to the ±X, ±Y, and 
±Z directions of the global coordinate frame, respec- 
tively. We have introduced the moment amplitude Mo = 
SJAiB(4a 2 +/3 2 — 2j 2 ) « 2.9fiB with the Lande factor gj — 
4/5. The sign x = ± represents the net pseudospin chi- 
rality of the tetrahedron, kt — \ ri r 2 r 3 ®t x 1 &r 2 x o> 3 , 
through the relation Ikt^/) = \fic\x^x,x' ■ Here, 
the summation over the sites r\,r2,r% on the tetrahe- 
dron T is taken so as they appear counterclockwise about 
the outward normal to the plane spanned by the three 
sites (Fig. [3] (b)). This kt gives the solid angle sub- 
tended by the four pseudospins (Fig.[3](c)). Note that the 
"2-in, 2-out" singlet state with the A\ g symmetry [26j . 
T, t =±(\ tX ) + \ tY ) + \ tZ ))/VG, is located at a high en- 
ergy AE ~ 7 J. The triply degenerate first excited states 
consist of only "3-in, 1-out" and "1-in, 3-out" , and are 
located at AE ~ J. In fact, any eigenstate of the single- 
tetrahedron Hamiltonian is described by either "3-in, 1- 
out" and "1-in, 3-out" configurations or "2-in, 2-out" 
and "4-in" / "4-out" configurations. Therefore, quantum 
effects of creating "3-in, 1-out" and "1-in, 3-out" con- 
figurations from the "2-in, 2-out" can not be taken into 



account in the single-tetrahedron analysis. 

To overcome this drawback, we numerically solve the 
model for the 16-site (N = 16) cubic cluster with the 
periodic boundary condition (Fig. Q] (c)). It is found 
that the ground states have a six-fold degeneracy labeled 
by the inversion (/) even(+)/odd(-) and the wavevector 
kx = (^,0,0), k Y = (0,2s,0), or k z = (0,0,^) with 
the energy ~ — 8.825J per tetrahedron. The states as- 
sociated with ki have a cooperative quadrupole moment 
defined on each tetrahedron, (Qt) — 0.0387Mq, where 



M^Sij with i,j 



X, Y. Z. Namely, the 



net magnetic moment Mt in each tetrahedron T points, 
for instance, to the ±Z directions with a higher proba- 
bility than to the ±X and ±Y (Fig. 13(d)). Such ferro- 
quadrupole order spontaneously breaking the three-fold 
rotational invariance about the [111] axes can occur in 
the thermodynamic limit. This ferroquadrupole moment 
Q T linearly couples to a lattice vibration: the four fer- 
romagnetic bonds and the two antiferromagnetic bonds 
should be shortened and expanded, respectively, leading 
to a crystal symmetry lowering from cubic to tetrago- 
nal accompanied by a compression in the direction of the 
ferroquadrupole moment (Fig. [3] (d)). This state shows 
both axial alignments of magnetic dipoles and a broken 
translational symmetry, and can then be classified into a 
magnetic an alog of a smectic (or crystalline) phase of liq- 
uid crystals [17] . Such magnetic quadrupole correlations 
in Pr2TA/207 could be probed by NMR experiments. 
Next, we calculate the magnetic dipole correlation, 

averaged over the six-fold degenerate ground states. This 
quantity is relevant to the neutron-scattering intensity 
integrated over the low-energy region below the crystal- 
field excitations from the atomic ground doublet Eq. ([TJ . 
Note that neutron spins do not couple to which rep- 
resents the atomic quadrupole. Figure 0] (a) shows the 
profile of S(q) for q = —(hhl). It exhibits intense peaks 
at (001) and (003), weaker peaks at (||0), and the min- 
imum at (000), as in the dipolar spin ice [8j, though the 
peaks are broadened by the quantum fluctuations in this 
ferroquadrupolar state. Besides, the nonzero q term in 
the Hamiltonian partially violates the ice rule and elim- 
inates the pinch-point singularity 
(111) and (002) in the spin ice [29 



which should be ex- 
amined by large system-size calculations. Note that our 
magnetic profile reproduces powder neutron-scattering 
results on P^S^Oy lfj that reveal the enhanced low- 
energy short-ranged intensity at q ~ — ~ 0.5 A -1 with 
a shoulder at q ~ — (Fig. BJb)). This experiment 
also shows a saturation of the quasielastic peak width 
~ 0.1 meV ~ J at 0.2 K [Hj]. Such large spin relaxation 
rate can be attributed to the appreciable quantum na- 
ture; large 5 and q in Eq. (2). These agreements support 
our scenario of the quantum melting of a spin ice. 

Now we concentrate on the ground states having the 



2JJ, [28| observed at 
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FIG. 4: (Color) (a) S{q)/Mo constructed from the local, 
nearest-neighbor, and second-neighbor correlations, (b) The 
theoretical curve (red) I(q)F(q) 2 with the form factor F(q) 
and the experimental data (green) of the powder neutron- 
scattering intensity on Pr2Sn207 at 1.4 K [lfj. 7(g) (blue 
curve) is the angle average of S(q)/Mo. (c) The magneti- 
zations (left) and energies (right) per site for the J-odd(-) 
ground state and the 7-even(+) state under the magnetic field 
H || (111), (d-f) Quadrupole correlations (QtQ",) between 
the tetrahedrons T and T' displaced by r = rx (d), ry (e), 
and rz (f) in the cooperative ferroquadrupolar state with 
the wavevector k z and (Qf z ) / 0. The matrix (QtQt ) 
in i, j has been diagonalized to yield two orthogonal forms 
of quadrupoles, Q rfJ , = -^mQt (M = 1j2). The shape of 
Q r i showing the dominant correlation amplitude is shown. 
Red/blue regions represent positive/negative values of Q r i. 



quadrupole moment > and the associated 

wavevector kz- The magnetic susceptibility is finite 
in the ferroquadrupolar state, as seen from the slope 
of the magnetization curve along the (111) direction in 
Fig. S] (c). This indicates a neg ative Tew as found in 
Pr 2 Zr 2 07 [l2j and Pr 2 Ir 2 07 14 1, and the absence of an 
internal magnetic field as in Pr 2 Ir 2 C>7 [l6j]- The magnetic 
field lifts the degeneracy due to the 7 symmetry. The 
ground-state (7-odd) magnetization shows a small step 
or dip around /is 77/ 'J ~ 1.5, in comparison with that of 
the 7-even excited state. This indicates that the structure 
develops upon cooling. These agree with the experimen- 
tal observation on Pr 2 Ir 2 07; M ~ at the metam- 
agnetic transition /zb77 c /J • 



14]. 



1.3 with J - 1.4 K 
Finally we spatially resolve the multipolar correlations 
within the cubic unit cell. Figures 0] (d), (e), and (f) 
represent quadrupole correlations {Q l ^Q 3 ^,) between the 
tetrahedrons T and T' displaced by rx = (0, §, § ), ry = 
(f,0, §), and r z = (|,§,0), respectively. There exist 
dominant ferroquadrupolar correlations shown in Figs. 0] 
(d) and (e), both of which favor ferroquadrupole mo- 
ments along the Z direction. They prevail over subdom- 
inant antiferroquadrupole correlations shown in Fig. [4] 
(f), and are responsible for the ferroquadrupole order 
(Qt Z ) 7^ 0- On the other hand, the chirality correla- 
tion (ktkt') is weakly ferrochiral between the tetrahe- 
drons shown in Figs. 0] (d) and (e), while it is strongly 
antiferrochiral between those shown in Fig. @] (f). This 



points to a geometrical frustration suppressing the chi- 
rality LRO in each fee sublattice of the diamond lattice 
formed by the tetrahedrons. Further studies are required 
for examining a possibility of a chiral spin liquid [3[ . The 
broken time-reversal symmetry without magnetic dipole 



LRO, reported in Pr 2 Ir 2 07 14], might be detected even 



in insulating magnets such as Pr 2 Zr 2 C>7 and Pr 2 Sn 2 C>7 
through magneto-optical Kerr-effect measurements. 

The proposed scenario of the quantum melting of 
the spin ice explains magnetic properties observed in 
Pr 2 TM 2 07. Effects of coupling of localized /-electrons to 
conduction electrons on the transport properties are left 
for a future study. The orbital motion of conduction elec- 
trons can flip the pseudospin-1/2. This could be an origin 
of the resistivity minimum observed in Pr 2 Ir 2 C>7 (l3| . 
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